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ABSTRACT 

Because of their large numbers, low mass stars may be the most abun- 
dant planet hosts in our Galaxy. Furthermore, terrestrial planets in the hab- 
itable zones (HZs) around M-dwarfs can potentially be characterized in the 
near future and hence may be the first such planets to be studied. Recently 



Dressing fc Charbonneaul ( 120131 ) used Kepler data and calculated the frequency 
of terrestrial planets in the HZ of cool stars to be 0.15l ; 06 P er s ^ ar f° r Earth- 



size planets 



.5 — 1-4 R©). However, this estimate was derived using the 



Kasting et al. 



( 119931 1 HZ limits, which were not valid for stars with effective 



tempera tures lower than 3 7 00 K. Here we update their result using new HZ lim- 



its from Kopparapu et al.l ( 120131 1 for stars with effective temperatures between 
2600 K and 7200 K, which includes the cool M stars in the Kepler target list. The 
new habitable zone boundaries increase the number of planet candidates in the 
habitable zone. Assuming Earth-size planets as 0.5 — 1.4 R®, when we reanalyze 
their results, we obtain a terrestrial planet frequency of 0. 48^24 an d 0.53l^?7 
planets per M-dwarf star for conservative and optimistic limits of the HZ bound- 
aries, respectively. Assuming Earth-size planets as 0.5 — 2 R e , the frequency 
increases to 0.511q2o P er s ^ ar f° r ^ ne conservative estimate and to 0.6llg 15 P er 
star for the optimistic estimate. Within uncertainties, our optimistic estimates 
are in agreement with a sim ilar optimistic estim ate from the radial velocity sur- 



vey of M-dwarfs (0.41±g;f|, iBonfils et al.l (1201 lh ). So, the potential for finding 



Earth-like planets around M stars may be higher than previously reported. 
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Introduction 



Advances in detection techniques and recent exoplanet s urveys have discovered terres- 



trial planets in the habitable zone (HZ) of their parent star (lUdrv et al. 



2010 



2012 



planets m tne nabitabie zone (t 
IPepe et al.lhpllai iBorucki et al.l 



2011 



Bonfi 



Tuomi et al.ll2012bl ; lAnglada-Escude et al. 



2007 



s et al. 2011; Borucki et al. 



Vogt et al. 



2012 



Vogt et al. 



2012) and it is expected that this number 



will greatly increase as time passes (IBatalha et al.ll2012l ). Classically, the HZ is defined as 



the circumstellar region in which a terrestrial mass planet (0.1 < M < 10 Mm), with 



able atmospheric conditions, can sustain liquid water on its surface (iHuang 



Kasting et al. 



1993 



Underwood et al.ll2003l ; ISelsis et al.ll2007bl ; iKaltenegger fc Sasselov 



1959 



Hart 



avor- 



1978 



2011 



Kopparapu et al.ll2013l ). Currently, more than 800 extra-solar planetary systems have been 
detected.], an d ~ 3000 additional candidate systems from the Kepler mission are waiting to 
be confirmed (IBatalha et al.ll2012f ). 



One of the primary goals of Kepler mission is to determine r? m , the fr equency of Earth- 
size planets in and near the HZ of solar-type stars (IBorucki et al.l 120 111 ). S uch estimates 
have b een made for pote ntially rocky planets using both radial- velocity (RV. iBonfils et al. 



(120111) ) and Kepler data (ICatanzarite &: Shad l201ll ; lTraubll2012l ; Dressing &: Charbonneau 
20131 ) . Fu rthermore, many other studies have estimated in general the terrestrial p lanet 



frequency (Howard et al.ll2010l ; Mavor et al.ll201ll ; iHoward et al.ll2012t ISwift et al.ll2012l ) and 
the consensus from these studies is that there are more low mass/smaller radii planets than 
high mass /large r radii ones. Moreo ver, the frequency increases towards lower radii, larger 
orbital periods (IHoward et al.ll2012l ). The estimates of the occurrence rate s of terrestrial 
planets with periods < 50 d ays ranged from 0.23ln'l 6 ( Howard et al.ll2010l ) around GKM 



stars using RV, to 0. 511° 05 (IDressing fc Charbonneau! 120131 ) around M-stars using transit 
detection. 

Specifically, to estimate 77©, one needs to know the boundari es of the HZs. Some stud- 
i es did estimate this numb er for M-dwarfs using transit (0 .15ln'rj v Dressing k, Charbonneau 



(120131 )) and RV (0.41±°;^, IBonfils et al.l (l201ll )) surveys. iTraubl (120121 ) estimated rj e to be 



0.34 ± 0.14 for Kepler FGK stars, but this is based on an extrapolation of data for planets 
with orbital periods shorter than 42 days. These estimates of the occurrence rates were 



based on 1-D radiative-convecti ve, cloud-free climate model calculat i ons by iKasting et al. 
Jl993l). Seve r al oth er studies (jUnderwood et al.l 120031 ; ISelsis et al.l l2007bl ) parametrized 



Kasting et al.l (119931 ) results to estimate relationships between HZ boundaries and stellar 
parameters for stars of different spectral types. Moreover, no uniform criterion of HZ limits 
were followed in applying them to calculate the occurrence rates, which can lead to com- 
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paring quantities that are not similar. For example, for the outer edge of the HZ, some 
studies ([Dressing fc Charbonneaul 120131 ) used the first CO2 condensation limit and others 
f lBonfils et al.ll201ll ) used Early Mars criterion. 



Recently iKopparapu et al.1 (120131) obta ined new, improved estimates of the boundaries 
of the HZ by updating Kasting et al.l ( 119931 ) model with new H2O and CO2 absorption coef- 
ficien ts from updated line- by-line (LBL) databa ses such as HITRAN 2008 (IRothman et al. 
20091) and HITEMP 2010 jRothman et alJboioh . They estimated that, for our Sun, a con- 
servative estimate of the inner HZ (IHZ) is at 0.99 AU and the outer HZ (OHZ) is at 1.70 
AU. These values represent the "water loss" (moist greenhouse) limit at the inner edge and 
"maximum greenhouse" limit at the outer edge. 



Kopparapu et al.l (120131 ) have also estimated HZ boundaries around other stellar spectral 
types, including M s tars, which are primary targets for ongoing surveys such as Habitable 
Zone Planet Finder (IMahadevan et al.l 120121 HPF) and MEarth (INutzman &: Charbonneau 
20081 ) to discover potential habitable planets. Furthermore, M-dwarfs are also increasingly 
becoming important as Kepler targets to find terrestrial size planets: The planetary orbital 
periods around these stars are shorter, the transit signal is larger, and there is a greater 
chance of detecting these planets than there is around a Sun-like star. Although M-dwarfs 
are fainter, the HZs are also closer for M-dwarfs, so it is likely that potential habitable 
planets may be discovered around nearby M-dwarfs in the near future with space-based 
characterization missions. 



In this study, we apply IKopparapu et al.l (120131 ) HZ limits to estimate the occurrence 
rate of terrestria l size (0.5 — 1.4 R m and .5 — 2 R ffl ) planets in the HZ of M stars. We base 



our analysis on Dressing fc Charbonneaul (120131 ) who updated stellar parameters of 3897 
low mass Kepler targ et stars with temperatures below 400 K using Darthmouth stellar 



evolutionary models (IDotter et al. 



2008 



Feiden et al 



20 111). T he outline of the paper is 



as follows: I n tf2\ we will revise Dressing fc Charbonneaul (120131 ) calculations with new HZ 
estimates of IKopparapu et al.l (120131 ) . discuss the implications in $3] and conclude in $U 



2. Estimate of habitable zone planet occurrence rate around M-dwarfs 



From lHoward et al.l (120121 ) and Dressing &: Charbonneaul (120131 ). the planet occurrence 
rate over a given period (P) and radius (R p ) range is given by: 



f(R P ,P) 



N P (R P ,P) 

E 



8=1 
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where dj is the semi-major axis of planet i, is the host star's radius of planet i, is 
the number of stars around which planet i could have been detected and N p (R p , P) is the 
number of planets with the radius R p and period P. The ratio ai/R*^ is the inverse of the 
probability of transit orientation, which is considered to take non-transiting geometries into 
the estimation of occurrence rate. 



In Table 2 of Dressing fc Charbonneaul (120131 ). the authors provide stellar and planetary 
parameters of ca ndidate KOIs that are con s idered to be in or near the HZ. To estimate the 



occurrence rate, 



Dressing fc Charbonneaul (120131 ) adopt the most conservative estimate of 



HZ limits from iKasting et al. f|l993h results: the 'moist greenhouse for the inner HZ (0.95 



AU) and the beginning of the CO2 condensation for the outer HZ (1.37 AU for the Sun). 
Furthermore, they consider 'Earth-size' as planets that are in the radius range 0.5 — 1.4 R e . 
Based on these definitions, iDressing fc Charbonneaul (120131 ) consider two KOIs to be in the 
HZ from their Table 2: KOI 2626.01 and 1422.02. They then calculate occurrence rate of 
Earth-size planets in the HZ as 0.15l o ; O g planets per star. 

This estimate may need to be updated, however, as the esti mated width and the po si- 
tion of the HZ has changed recently, following new calculation by Kopparapu et al. ( 2013 )^1. 
According to these authors, HZs are farther out from their star than what has been cal- 
culated previously. This will h ave a significant effect on the occurrence r ate derived by 



Dressing Sz Charbonneaul (120131 ). Also, The "1st CO2 condensation" limit of IKasting et al. 



(119931 ). should n ow be disregarded, as it has be en shown that CO? clouds generally warm a 
planet's climate (IForget fc Pierrehumbertlll997l )). Here we use the Kopparapu et al.l (120131 ) 
HZ limits to derive revised estimates of occurrence rates of potentially habitable planets 
around low- mass stars. In order to put a lower and upper bound on 77®, we calculate two 
values for t he occurrence rate based on the conservative and optimistic estimate of HZ limits 



as given m 



Kopparapu et al.l (120131 ): 



(1) In a conservative estimate, the inner edge of the HZ is determined by the "moist- 
greenhouse" limit which is derived by assuming a fully saturated troposphere and negligible 
cloud feedback. Neither assumption is likely true in reality, but it is difficult to improve on 
this with a 1-D climate model because such models cannot accurately simulate clouds or 
relative humidity. The outer edge of the HZ is determined by the "maximum greenhouse" 
limit where a CO2 dominated atmosphere can produce maximum amount of greenhouse 
warming. Here also, the radiative warming by CO2 clouds is neglected hence the limit is a 
conservative estimate. 

(2) In an optimistic scenario, the inner edge of the HZ boundary can be obtained by the 



^http : //www3 . geosc . psu . edu/~rukl5/planets/ 
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"recent Venus" limit which is based on the observations of Venus by Magellan spacecraft, 
suggesting that liquid w ater has been absent from the surface of Venus for at least 1 Gyr 
( jSolomon fc Headlll99ll ) or earlier. The Sun at that time was ~ 92% of the present day lumi- 



nosit y, according to standard stellar evolutionary models (IBaraffe et al.lll998t iBahcall et al. 



2001 



See Table 2). The current solar flux at Venus distance is 1.92 times that of Earth. 
Therefore, the solar flux received by Venus at that time was 0.92 x 1.92 = 1.76 times that of 
Earth. This empirical estimate of the inner HZ edge in our Solar System corresponds to an 
orbital distance of d = (1/1.76) ' 5 = 0.75 AU for the present day. Note that this distance is 
greater than Venus' orbital distance of 0.72 AU because the constraint of surface water was 
imposed when the Sun was fainter. The outer edge optimistic estimate is the "Early Mars" 
limit based on the observation that early Mars was w arm enough for liquid water to flow 
on its surface (jPollack et al.l 119871 ; iBibring et al.ll2006l ). Assuming the dried up riverbeds 
and valley networks on martian surface are 3.8 Gyr old, the solar luminos ity at that time 
would have been ~ 75% of the present value (See Eq.(l) in iGoughl (I198ll ) and Table 2 in 



Bahcall et al.l (120011 )). The present-day solar flux at Mars distance is 0.43 times that of 
Earth. Therefore, the solar flux received by Mars at 3.8 Gyr was 0.75 x 0.43 = 0.32 times 
that of Earth. The corresponding empirical OHZ limit for our Solar System today, then, 
would be d = (1/0.32) ' 5 » 1.77 AU. 

The optimistic HZ limits can be extrapolated to other stellar types with effective tem- 
peratures between 2600 K to 7200 K, by scaling them wi th the corresponding values of 
conservative HZ limits, as shown in iKopparapu et al.l (120131 ) . 



To derive new rates from Eq.(pQ), we need two quantities: (1) ai/R+j for the planets 
in the HZs and (2) N* t i, the number of stars around which a planet that has the same size 
and receives the same insolation as those KOIs in the HZ could have been detected. The 
photo metrica lly derived values of aj/R + j are given in Table 5 of Dressing &: Charbonneau 



(120131 ). But Dressing fc Charbonneau! (120131 ) use the calculated ai/R+j determined from 
the period of the planet and mass of the star (Dressing, private communication). We use 



the N + j values provided within the period-radii cells of Fig. 15 of iDressing Sz Charbonneau 



(120131 ). which gives the number of stars around which a planet from the center of the grid cell 
would have been detected with a signal to noise ratio above 7.1 a. This should still give us 
nearly the same occurrence rate, or an underestimate of th e actual value (see next section) . 
For example, for the two KOIs (2626.01 and 1422.02) that bressing fc Charbonneau! (l2013h 
consider to be in th e HZ, when we use the center o f the grid cell iV*^ numbers (1822 and 
872) from Fig. 15 of IDressing fc Charbonneau! (120131 ) (instead of the iV*^ values provided in 
section 5.7) and use the period determined di/R^^ as they did, we were able to reproduce 
their value of 0.15. Thus, our new occurrence rate (which uses center of the grid cell numbers) 
determined here probably is close to a value from a more rigorous estimate. 
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We will first calculate t)q using IDressing fc Charbonneau ( 2013 ) radius range of 0.5 — 1.4 



R ay Assuming conservative HZ lim itajj from lKopparapu et al. 



(12013 1. four KOIs from Table 2 



of lDressing fc Charbonneaul ( 120131 ) should be in the HZ based on the insolation fluxes (FigfTJ) : 
K OI 1686.01, 2418.01, 2626 01 and 1422.02. We use the photometric a^R^ from Table 5 
of lDressing &: Charbonneaul (120131 ). rather than the calculated a^/i?*^ from the period and 
stellar mass. For a typical planet candidate, the photometric ai/R+j is 85% of the a*/ ' R*,i 
from the period and stellar mass. Therefore, our value tends to be at the lower end of the 
occurrence rate estimate. The corresponding aj/i?^ for these candidates from Table 5 of 
Dressing fc Charbonneaul (boia) are: 102.482, 116.837, 36.283 and 51.985, respectively. The 
corresponding N+j values from center of the grid cells of Fig. 15 are: 353, 994, 1822 and 872, 
respectively. Using Eq.flT]), we get the conservative estimate of 77© for low- mass stars to be 
0.48±[};24 per star. 

An optimistic estimate of the occurrence rate can als o be derived based on the recent 
Venus and Early Mars li mit from iKopparapu et al. (120131) results. Assuming these limits, 
six KOIs from Table 2 of lDressing fc Charbonneau ~J2013|) should be in the HZ (FigflJ: KOI 
1686.01, 2418.01, 2626.01, 1422.02, 2650.01 and 886.03. The corresponding photometric 
determined Qj / R*^ for these candidates from Table 5 of IDressing &: Charbonneaul (120131 ) 
are: 102.482, 116.837, 36 283 51.985, 54.052 and 39.246. The values from Fig. 15 of 
Dressing fc Charbonneaul (120131 ) are: 353, 994, 1822, 872, 1822 and 2336. Using Eq.©, the 
optimistic estimate of the occurrence rate of Earth-size planets in the habitable zones around 
low-mass stars is 0.53lo 17 P er star. 

Instead of assuming 0.5 — 1.4 R© as 'Earth-size', we will also calculate occurrence rate 
extending the radius range from 0.5—2 R ffi . Planets with radius > 1.4 R ffi are thought to have 



(Seager et al. 


2007 


„ ^ — , 

Rogers et al. 


2011; 


Lopez et al. 


2012) 



add KOI 854.01 in the conservative rate estimate (the total number of KOIs in the HZ is 
then five) and KOI 250.04 in the optimistic rate estimate (the tot al HZ KOIs is eight). The 
corres ponding cii/R+j and for these additional candidates from IDressing Sz Charbonneau 



( 120131 ) Table 5 and Fig. 15 are: 90.045 and 2991, respectively, for KOI 854.01; 157.259 and 



3287, respectively, for KOI 250.04. 



Using this expanded definition of Earth-size (0.5 — 2 R©), a conservative estimate of the 
irrence rate of Earth-size planets in the HZs e 
optimistic estimate on the occurrence rate is 0.61" 



occurrence rate of Earth-size planets in the HZs around M-dwarfs is 0.5llg 20 P er star. An 



"q;?5 per star. 



3 ie., Moist-greenhouse limit at the inner edge and maximum greenhouse limit at the outer edge 
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Discussion 



The occurrence rate estimates derived in the previous section indicate that terrestrial 
size planets in HZs around low- mass stars are more frequent than previously thought. Couple 
of caveats are to be noted: 

(1) The KOIs are not confirmed planets. The calculat ed false positive ra tes for Ke- 
pler candidates in our rel evant radius bin range from 12.3% ( jFressin et al.l 120131 ) to ~ 10% 
( iMorton fc Johnson! |20 111 ). Therefore, we may be overestimating 77®. Furthermore, the un- 
certainties on the fluxes are large (Fig. [TJ), so some may not be i n the HZ (2) On the other 



Dressing fc Charbonneau 



Dressing fc Charbonneau 



hand. We use photometrically derived / R+j given in Table 5 of 
( 2013 ). This may underestimate our occurrence rate because 
( 120131 ) use calculated a^/i?*^ from orbital periods and stellar mass (Dressing 2013, private 
communication) and the photometric values are ~ 85% of the derived values. So the net 
effect from points 1 and 2 may not change significantly our estimate of 77®. 

(3) As mentioned above, we do not calculate number of stars around which a planet that 
has the same size and receives the same insolation as the ones in the HZ could have been 



detected; we instead use the cell numbers from Fig. 15 of IDressing fc Charbonneau! ( 120131 ). 
which gives the number of stars around which a hypothetical planet that has the same radius 
and period as the center of the grid cell would have been detected with a signal to noise ratio 
above 7.1cx. There is no reason to assume that the KOIs considered here to be in the HZ are 
at the center of the grid cell (in fact, they are not). This offset probably overestimates our 
value of N* j in Eq. (Tj[|) . 



For example, in section 5.7 of IDressing fc Charbonneau! (120 131 ). they calculate N+i values 
of 2853, 813 and 2131 for KOIs 854.01, 1422.02 and 2626.01, respectively. These numbers 
are generally lower in value than the center of the grid cells values we use from Fig. 15 
(2991,872, 1822, respectively), except for KOI 2626.01 (an increase of 17%). If this is the 
general trend, i.e, if we are overestimating systematically, then our calculated occurrence 
rates of terrestrial planets in the HZ around M-dwarfs can be considered as a lower bound 
to the actual value. Nevertheless, we tested to see how much the rate would change if we 
change iV*^. It should be noted that N* t i (at same insolation) is most likely smaller than JV*^ 
(at same period) for the planet candidates in the HZ because most of those planet candidates 
orbit stars that are cooler than the typical Kepler M dwarf. Since there is one KOI (2626.01) 
for which we underestimate by ~ 17%, and could potentially decrease our occurrence 
rate estimate, we rounded off and added 20% of the respective N^i values to all the KOIs 
that are in the HZ. After performing this calculation, for the radius range 0.5 — 1.4 R®, the 
conservative estimate changed from 0.48^024 to 0. 41^20 P er star. The optimistic estimate 
changed from 0.53+oi7 to 0.44toiI per star. In the extended radius range for Earth-size 



the conservative estimate changed from 0.511° 20 per star to 0.431q 175 and 



(0.5 - 2.0 R e „ ^ — , „ c — -0.20 

the optimistic estimate changed from 0.6llo!?5 P er star to O^Iq^. 



Note that even though adding 20% to iV*^ lowers our occurrence rate, we are underes- 
timating ai/ R+^i on an average of 15% or mo re. Furthermore, we showed that we were able 
to reproduce Dressing fc Charbonneaul (120131 ) estimate even if we use center of the grid cell 
v alue for N+ j. Therefore , a mo re rigorous analysis to determine the occurrence rate from 



Dressing fc Charbonneaul (120131 ) would likely produce a similar result as our's. 



Bonfils et all ( 1201 ll ) studied 102 southern nearby M dwarfs using ESO/HARPS spec- 
trograph and obtained the frequency of terrestrial mass planets (1 — 10 M ffi ) in the HZ to 
be 0.411q i3- Note that they consider two planets to be in the HZ, Gl 581d and G l 667Cc , 
using optimistic limits of the HZ (recent Venus and Early Mars) from lSelsis et al.l ( 12007bl ). 
If we use also use optimistic HZ estimates, we find that both the planets are in the HZ. 
Our optimistic es timate of r anges from 0.53tQ;57 to 0.6llo!i5 P er star, depending on if 
the mass range in lBonfils et al.l ( 1201 ll ) is applicable to terrestrial sizes of 0.5 — 1.4 R® radius 
or 0.5 — 2 R® radius. Nevertheless, our estimates of the occurrence rates of HZ terrestrial 
planets around M-dwarfs are in good agreement with independently derived RV estimates, 
within the error bars. 



Conclusions 



The purpose of our analysis is to outline the signi ficant increase in the terre s trial p lanet 



occurrence rate in the HZs of M-dwarfs compared to 



Dressing fc Charbonneaul (120131 ). We 



have obtained revised estimates based on iDressing fc Charbonneaul (120131 ) estimates of 3987 
Kepler M-dwarfs t hat a re cooler than 4000 K in Q1-Q6 data. Applying new HZ results from 
Kopparapu et al.l (120131 ) to Earth-size planets (0.5 — 1.4 R©), we calculate that a conservative 



estimate of the occurrence rate of Earth-size planets in the HZs around M-dwarfs is 0.48lo24 
planets per star. The optimistic estimates indicate that the occurrence rate increases to 
0.53jlo 17 planets per star. If we extend the definition of Earth-size to planets in the radius 
range 0.5—2.0 R®, the conservative estimate increases to 0.51 1^20 P er s ^ ar ' ano ^ the optimistic 
estimate increases to 0.6lio!i5 P er star. As discussed in our paper, our results probably are 
close or underestimate the actual occurrence rate of Earth-size planets in the HZs around 
M-dwarfs. Furthermore, our optimistic value of 77® qu oted above is i n agr eement (within 



uncertainties) with a similar estimate, 0.411q; 13 , from iBonfils et al.l ( 20111 ) ESO/HARPS 



survey of 102 M-dwarf stars indicating that the frequency of terrestrial planets in the HZs 
of M-dwarfs may be higher than previously reported. 
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Fig. 1. — Incident stellar flux on a planet as a function of stellar effective temperature, T e ff. 
The green shaded region is the conservative HZ. The optimistic HZ limits are recent Venus 
(solid red curve) and Early Mars (solid blue curve). Tw o of th e terrestrial-size KOIs (1422.02 
and 2626.01) that are in the iDressing fc Char bonneaul (120131 ) HZ are also shown. 
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